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Abstract 

We provide detailed analysis of strong field ionization of degenerate valence p orbitals by circu- 
larly polarized fields. Our analytical approach is conceptually equivalent to the Perelomov, Popov, 
and Terent'ev (PPT) theory and is virtually exact for short range potentials. After benchmarking 
our results against the PPT theory for s orbitals, we obtain the results for p orbitals. We also show 
that, as long as the dipole approximation is valid, both the PPT method and our results are gauge 
invariant, in contrast with widely used strong field approximation (SFA). Our main result, which 
has already been briefly outlined in [I. Barth and O. Smirnova, Phys. Rev. A 84, 063415 (2011)], 
is that strong field ionization preferentially removes electrons counter-rotating to the circularly 
polarized laser field. The result is illustrated using the example of Kr atom. Strong, up to one 
order of magnitude, sensitivity of strong field ionization to the sense of electron rotation in the 
initial state is one of the key signatures of non-adiabatic regime of strong field ionization. 

PACS numbers: 42.50.Hz, 32.80.Rm, 33.80.Wz 



I. INTRODUCTION 



The analysis of ionization in strong low-frequency laser fields is often based on adiabatic 
approximation. In this approximation, ionization is treated as tunneling through a static (or 
quasi-static) barrier created by the binding potential and the voltage drop due to the electric 
field of the laser pulse. This picture implies that the electron does not see the oscillations 
of the low-frequency laser field during ionization, i.e. tunneling happens "faster" than the 
oscillation of linear or the rotation of circular field. Formally, this picture corresponds to 
the limit 7 < 1, where 7 = 2I P u/£ is the Keldysh parameter [1], I p is the ionization 
potential, u is the laser frequency, and £ is the strength of the laser field. 

However, for typical experimental conditions, both for linear and for circularly polarized 
fields [2HZ], the Keldysh parameter is often in the non-adiabatic tunneling regime [El E], 
i.e. 7 ~ 1. Therefore, the adiabatic-based interpretation of these experiments is questionable. 
In particular, we have shown in Ref. pU] that for strong field ionization in circularly polarized 
laser fields, the sense of electron rotation becomes significant already for 7 < 1, i.e. even 
when using longer-wavelength laser radiation than the standard 800 nm, e.g. 1300 nm as in 
Ref. [5]. As a consequence of non-adiabatic effects, the counter-rotating electron can have 
up to one order of magnitude larger ionization rate than co-rotating, depending on the laser 
field parameters. Our theoretical prediction has now been confirmed by the experiment [TT] . 

The goal of this paper is to expose our calculations and discuss the approximations we 
have used in deriving simple formulas for ionization rates presented in Ref. [10J. We follow 
the theory developed by Perelomov, Popov, and Terent'ev (PPT) [T2l [13] for short range 
potentials and apply it to p orbitals. Including effects of the long-range Coulomb potential in 
a standard way [T3J [15] and including non-adiabatic Coulomb effects [HI [T7] do not change 
our conclusions. The Stark shift of the initial state is not included in our analysis, but it can 
be calculated separately [TJ HE] and used to correct the field-free ionization potential used 
in the present calculation. 

The key advantage of the PPT approach is its gauge invariance, which is discussed below 
in the paper. The violation [T9H2T] of gauge invariance in strong field approximation (SFA) 
leads to both technical and conceptual problems [22]. In particular, for the ionization of 
p orbitals by strong circularly polarized fields, the SFA yields inconsistent results in both 
gauges [231 121] , which contradict experimental measurements of either ionization yields [11] 
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or electron spectra (see e.g. Refs. [2H1J[6J[7]). Note that in linearly polarized fields the length 
gauge SFA and the PPT theory yield equivalent results for short range potentials [2S]. The 
results of length gauge SFA and the PPT are different if long range effects are taken into 
account, with length gauge SFA leading to incorrect prefactor of ionization rate pp. The 
deficiencies of velocity gauge SFA are well documented [T9H2T] and are significant even for 
short range potentials, e.g. velocity gauge SFA predicts identical total ionization rates from 
p + and p_ orbitals [23J. 

Our analysis reveals that optimal quantum trajectory, which minimizes electron action 
under the barrier, corresponds to initial non-zero lateral velocity in direction opposite to the 
rotation of the laser field. The weight of this trajectory is determined by the direction of 
electron current in the initial orbital and is higher for p_ orbitals in case of right circular 
polarization of the laser field. 

Finally, we show that non-adiabatic dynamics of strong field ionization leads to non-trivial 
rotational dynamics of the hole left in the ion. This dynamics leads to the generation of elec- 
tronic ring currents in ions [26, 27J, and the coherence of this dynamics can be probed with 
attosecond time-resolution using attosecond transient absorbtion technique demonstrated 
recently in Ref. [28J. 

This paper is organized as follows: In Section II, we derive analytical formulas for the 
ionization rates based on the PPT theory. We first benchmark our results for s orbitals 
against the PPT results [12] . We then derive the results for po and p± orbitals. In Section 
III, we apply these formulas to the strong field ionization of the Kr atom. Section IV 
concludes this work. 

II. THEORY 

A. Ionization model in circularly polarized laser fields 

The PPT formulas for the atomic ionization rates were derived for the strong field ioniza- 
tion in linearly, circularly, and elliptically polarized laser fields [T2| IT3] . However, for circular 
and elliptical polarizations, there are formulas only for s orbitals. In this section, we derive 
the analytical formula for the ionizaton rates in circularly polarized laser fields also for p m 
orbitals with azimuthal quantum numbers m = 0, ±1. The right (+) or left (— ) circularly 
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polarized laser field is defined as 



E±(£) = S (cos(u)t) e x ± sin(W) e y ) , (1) 

which is connected with the vector potential 

A±(i) = — A (sin(atf) e x =p cos(c<;t) e^) (2) 

by the relation E±(i) = —dA±(t)/dt, where £ is the electric field amplitude, Aq = £/u> is the 
velocity amplitude of the electron oscillations in the laser field, and u is the laser frequency. 

We assume that the electron ionization from the valence p shell is described by the 
time-dependent Schrodinger equation (TDSE) in single active electron model within dipole 
approximation 



• 9 , , , 



v 
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-£ + y(r)+r.E±(f) 



V±M), (3) 



where V(r) is the effective potential and the atomic units are used throughout in this work. 
The exact solution of this TDSE is the integral equation for the time-dependent wavefunction 
ip±(r,t) starting at time t 



V±(r,t)= / dT f G ± (T,t,T f ,t )iP ± (T f ,t Q )-i f dU [ dr f G ± (r,ty,U)V(T f )^ ± (T f ,t i ), (4) 

J J t J 

where 

G ± (T,tJ,U) = 9 4 JZ ^T- f dke^>-^) r V^ v± M 2 * (5) 
(2ny i J 

is the Green's function of the electron for motion in a circularly polarized field, 

v±(t)=k + A±(*) (6) 

is the instantaneous electron velocity, and k is the final momentum observed at the detector. 
Moreover, we divide k = ki + kj_ into two components kii = k x e x + k y e y and kj_ = 
k z e z , which are parallel (kn || A±(t)) and perpendicular (kj_ _L A±(t)) to the laser field, 
respectively. 

The first term of Eq. Q does not contribute to the ionzation rate, because it describes 
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the smearing out of the initial state [12] . As in PPT theory, the main approximation of this 
theory is the neglect of the disortion of the initial wavefunction ip±(r',ti) by Stark effect 
prior to ionization at time tj, i.e. we replace the exact wavefunction 4>±{ r 'iti) on the right 
side of Eq. Q by the wavefunction of the bound orbital for the free atom {pi m (r')e lIpti with 
quantum numbers I, m and ionization potential I p . Using the field-free TDSE, the term 
V{y')iP±{y' , U) is replaced by 

V(r')w m (r')e iIpti = ~ 2/ P )^ m (r')e^. (7) 

As already described in Ref. [T2] in detail, the difference between two wavefunctions ip±(r', ti) 
and ipi m (r')e lIpti is small for short-range potentials, i.e. the potential V(r) falls more rapidly 
than the effective Coulomb potential ~ 1/r. However, the Coulomb corrections can be 
introduced using standard recipes [T^l ITS] involving the time-integration of the Coulomb 
potential along the optimal trajectory. In this work, we use the short-range potential and 
will include Coulomb corrections in our future work. 



Furthermore, we assume that the laser field is turned on at to ~^ — 00 adiabatically. Then, 
using the momentum representation of the wavefunction 



^™( k ) = 7^W2 J dre (8) 

and the abbreviation 

<f> lm (v±(t)) = l -{v ± {t) 2 + 2I J ,)ft m (v±(*)), (9) 
we get the approximative electron wavefunction from Eq. (|4]) 

V>±(M) = TTT^R f dt ^ eUpU I dke^ r e-^^ ±iT)2dT <p lm ^±(U)). (10) 

For circularly polarized fields, it is advantageous to use cylindrical coordinates (p, 0, z) in- 
stead of Cartesian ones (x,y,z) in coordinate space related by x = pcos0 and y = psin0. 
Similarly, we introduce cylindrical coordinates (k p , 9, k z ) in momentum space with relations 
k x = k p cos 9 and k y = k p sin9, thus k 2 p = k 2 x + k 2 = k 2 and k 2 = k 2 + k 2 = k 2 + k\. With 



Eqs. (|2]) and ([6]), two exponents in Eq. (10) are expressed as 



iv±(t)r = if ± (k p ,9,(f),t)p + ik z z (11) 

and 



\ f v ± (r) 2 dr = - l -(k 2 + Al) (t - tl ) - ik(£t(f) - £fc(i0), (12) 



where 

f±{k p , 9, (ft, t) = k p cos(# — 0) — A sin(wt =F 0) (13) 

and 

Z±(t) = E ± (t)/u 2 . (14) 



B. Derivation of the formula for the time-averaged ionization rate 

We follow the derivation of the formula for the ionization rate in Refs. [T21 [T3] based on the 
PPT approach and repeat it here for clarification and only for the case of circular polarization 
(e = 1). The time-averaged ionization rate w±(£,cj) is equal to the time-averaged radial 
flux at the infinity p — > oo, i.e. 

w±(£,u>) = lim J±(p,t). (15) 

The radial flux J±(p,t) is evaluated as the integral of the radial component of the flux 
density j P ±(r, t) over a cylinder of radius p with its axis along the propagation z-axis of the 
circularly polarized laser field, i.e. 

/OO /*27T 
dz j d(f)jp±(p,(f),z,t), (16) 
oo J0 

where j p ±(r,t) is defined as 

jp±M) = \ U±(T,t)^r±{T,t)-r±(T,t)^M^t)\ . (17) 
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Inserting Eqs. (10)-(12) into Eq. (17), we get 



Jp±M) = 2^ J dk! y dk 2 e i ( k - kl )( r -« ± W)(/ ± (A; lp ^ 1 ,0,t) + / ± (A; 2p ^ 2 ,0,t)) (18) 
f dt li eK fc ?+^+ 2 ^)( t - t ")l^(k 1 ,t li ) /"* dt M e-K^+^+ 2i Sp)(*-^) F± ( k2) t 2 .) ) 

J ~oo J— oo 



where the function 



F ± (k J t) = lm (v ± (t))e ik ^W 



(19) 



contains terms with complicated, but periodic time-dependence. Expanding F±(k, t) into 
the Fourier series 



F±(k,t) = ^ F n± (k,o;)e- 



inujt 



(20) 



with the Fourier coefficients 



1 



F n± (k, tu) = — I d(ut) F ± (k, t)e inw ' 

Z7T 



and carrying out time-integrations in Eq. (18), according to (a G M, 5 > 0) 



dt'e ±ia ^ = lim 



(21) 



(22) 



s->6 a ± z5' 

yields the final formula for the radial component of the flux density (5 > 0) 
Jp± (r,t) = lim^-p J dk! y ^^^-^^(^(/i^^^^^ + Zi^,^,^^) (23) 



E F n 1± (ki,u; 

L=— oo 

oo 

^2 F n2 ±(k 2 ,u 



+ — + ip - niw + iS 



Tl2= — OO 



— + — 

2 2 



n -1 



2 + ^ + I p -n 2 u-i5 



-i(n2—ni)u}t 



This expression is then inserted into the formula for the the radial flux J±(p,t), Eq. (16). 
With 



i(k 2 - ki)r = ip((k 2x - k lx ) cos <p + (fc 2y - k ly ) sin 0) + i(k 2z - k lz )z, 



(24) 
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the ^-integration is easily carried out, i.e. 

/oo 
diZ e i(ka-ki)p = 2ir6(k 2z - k lz )e ipak2x - klx)cos * +{k2y ~ kl y )s ' m ^. (25) 
-oo 



But the evaluation of the (^-integration is challenging. Using Eq. (13), the Euler's formula, 
and the substitution |k 2 || — kiy | sin 0' = {k 2x — k lx ) cos <fi + {k 2y — k\ y ) sin 0, the ^-integral is 
evaluated as (see Appendix 1) 

p2tt 

/ rf eV ((fe,-^)cos0 + ( fc , y - fcl!/ )sin0) (/±(A;i/j5 g u 0; t) + /±( ^ ^ 0; t)) (26) 
Jo 

= 2m {k 2 - k 2 - 2k 2p A sm(ut =p 9 2 ) + 2k lp A sm(ut T Jl ^ k ^~ k ^\ 

l k 2|| — | 

where J n (x) is the Bessel function of the first kind, cf. Ref. [13]. Then, we carry out the 
simple integration over k\ z and get the result for the radial flux 

J±(p, t) = lim % - J dki|| J dk 2 \\ e-^ k 2|[- k ui)e±(*) (27) 

( fc 2 P - - 2k 2 P A sm(ut =F 6» 2 ) + 2k lp A sin(wt =F 0i)) 

/OO 00 
cifc, ^ K 1± ((ki||, k z ), oo) [k 2 lp + k 2 z + A 2 + 2I p -2 ni oo + id]' 1 
-oo _ _ „ 



pjl (p 


k 2 ||-ki|||) 




k 2 || - kiy 





F n2± ((k 2h k z ),u) [k\ + k 2 z +A 2 + 2I P - 2n 2 u - i5] 1 e-^-^' 

ri2=— oo 



To obtain the limit of the radial flux at the infinity p — > oo, we apply the relation for the 
arbitrary function g(k||) (see Appendix 2) 

hm J dk\\ <7(k„) = 2n J dk„ ff(k||)*(k||), (28) 

cf. Ref. [13], carry out the integration of the radial flux over ki y , and use Eqs. ( [Tip , and 
the substitutions /c-t = k 2x ± fci^. Then, the intermediate result is 



oo , , 

z An K 



lim J ± (p,t) = \mvi I dk_ 5{k^)k^e — ^ cos M) (29) 

p— >oo 5— >0 



— oo 

oo oo 



/OO /'OO - v - v /*oo 

dk y dk z J2 e~ i{n2 - ni)ujt dk + h±{k + ) 

m=-oon 2 =-oo J ~°° 
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(k+ - k-) 2 + ky + k; + A 2 + 2I P - 2n x u + id 
(k+ + k_ Y + k 2 y + k 2 z + Al + 2I P - 2n 2 u - i5 



where the analytical function h(k + ) is defined as 

h±{k+) = (k+ - 2A sin(wt)) F* l± ( ( k+ k ~ ,k y ,k z J ,u) F n2 ± ( ( k+ * k ~ , k y , k z ) ,u ) . 



(30) 



It is now shown in Eqs. (29) and (30), that the ionization rate depends on the sense (±) of 



circular polarization only in the function F± (k, t) , Eq. ( 19 ) . By further deep analysis, the k 



integrand in Eq. (29) has four poles whose locations on the complex plane and corresponding 
residues depend particularly on n\ and n 2 . For k- = and ^ n\ 7^ n 2 7^ 0, all four residues 



are finite, thus the fc_-integral in Eq. (29) would be zero due to the appearance of the factor 
k_ in the integrand. Therefore, the condition for the number of photons n = rii = n 2 must 
be satisfied. Furthermore, for 2nu < k 2 + k 2 + A\ + 2I p , there are only two residues that 
could contribute to the ionization rate, but in the limit k- = these residues are opposite. 
Therefore, we consider only the case 2nu > k 2 + k 2 + Aq + 2I P , where only two residues 
contribute to the ionization rate. By the way, we denote the quanity no as the minimal 
number of photons required for ionization in circularly polarized laser fields, i.e. 



n 



A 2 
2cu 



UJ 



2U V + L 



UJ 



(31) 



where U p = Aq/4 is the pondermotive potential. Comparing to the case for linearly polarized 
laser fields, the mean kinetic energy of the electron in a circularly polarized laser field 
Aq/2 = 2U P is twice as much. Applying the residue method for the £; + -integral and evaluating 
the /c_-integral, we get the expression for the radial flux at the infinity (see Appendix 3) 



lim J±(p,t) = 7r 



p— >oo 



where 



n>no 



elk*. 



(lkz ^ (Vg - g -kp ~ h ± {-2^kj - kj - g) ^ (32) 



Z^»2 Z^«2 Z^"2 

K n H y K z 



k 2 

2 



{n — no)uj. 



(33) 
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Using Eq. (30) for = and n = n\ = n 2 , i.e. 

h±(k + ) = (k + - 2A Q s\\i(ujt)) 



F n ± ( ( — ,k y ,k z ) ,u 



(34) 



time-averaging over a laser cycle, and using the relations for the 5 function 5(a 2 — x 2 ) = 
[S(a — x) + 5(a + x)]/(2\a\) and 6 (ax) = 5(x)/\a\, we obtain the final formula for the 



ionization rate from Eq. (15) as a sum over multiphoton channels 



and 



w±(£,u) = ^2 w n± (£,uj), 

n>no 



w n ±(£,u>) = 2tt / dk5 [ y - S ) \ Fn± (k, W )' 2 



(35) 



(36) 



which coincide exactly with Eqs. (13) and (14) of Ref. [13] for circular polarization, respec- 



tively. In Eq. (36), we recognize that there is the conservation law, namely k = k n , that 



means that the electron kinetic energy is equal the photon energy minus the mean electron 



energy in a circularly polarized laser field and the ionization energy, cf. Eqs. (31) and (33), 
i.e. 



k 2 k 2 

Y = f = nu-2U P -I p . 



(37) 



C. Gauge invar iance 



The length gauge was used in the derivation as in the original PPT approach. However, 
we note that the result for ionization rate is independent of the gauge. Rewriting Eqs. (J3|- 
(|5| using the velocity gauge yields substitution of the original wavefunction ip±(r,t) in the 



length gauge by ip±(r,t)e * A ±(') r in Eq. (10). Thus, in right hand side of Eq. (10) the term 



e iv ± (t)r - g ^j ien replaced by e iV± ^ r e jA ±(*) r = e lkr and the function in Eq. (13) is therefore 
time- independent, i.e. f±(k p , 9, 



k p cos( 



). Following the derivation in the previous 



section, Eqs. (26) and (30) do not have any time-dependent terms anymore, yielding the 



time- independent radial flux in Eq. (32). Therefore, time-averaging over a laser cycle is 



unnecessary in the velocity gauge, yielding the same result for the ionization rate as in 



10 



Eqs. (35) and (36). 



D. Derivation of the formula for the probability of the n-photon process 



The function \F n ± (k, u)\ 2 in Eq. (36) describes the probability of the n-photon process 
in circularly polarized fields, which is derived in this work not only for s- but also for all 



atomic orbitals, thus beyond the derivations in Refs. [121 03] • With Eqs. (19) and (21), we 
start with the general formula for the probability of the n-photon process dbt k — k n 



\F n ±(k,u)\ k=kn 
where the action S±(k, t) at k = k n is 



CO 

4^2 



tt/lu 



dt <j>im{v±{t))e 



— ir/u> 



(38) 



S± (k, t) | k=kn = k (£* (f ) - k (0) ) | k , + nut. 



It is also obtained by the well-known expression for the action 



(39) 



i r* 

S±(M) = -J v±(r) 2 dr + I p t, 



(40) 



the conservation law (37), and the relation 
1 « 



v ± (r) 2 dr 



k 2 



2U p )t + k($ ± (t)-Z±(0)), 



(41) 



cf. Eq. (12). Using the saddle point method applied for to <ti I p and to U p , we obtain the 



simple expression for the probability of the n-photon process 

<Mv ± (t,)K 5±(kA) 



to 



\F n± (k,u)\ 2 k=kn = — 



2tt 



(42) 



k — kn 



It means that the integral in Eq. (38) is accumulated mostly in the small region around the 



so-called (complex) ionization time £j which is uniquely linked to k and determined by the 
saddle point equation 

v±(t,) 2 



d_ 
Of 



S±(k,t) 



k — kft Aj — ij 



+ nuo = 0. 



(43) 



k — kn ,t — t i 



n 



With 



(44) 



the saddle point equation (43) is rewritten as 



sm(uti =F 6) = Xn{k z ), 



(45) 



where 



Xn(K 



nuo 



nuj 

k 2 x ~ A^K 



> 



Xn{k z = 0) =: Xr, 



Using Eqs. (31) and (33), the variable Xn is reexpressed as 



n 2 oj 



Xr. 



n 2 (l + 7 2 ) 



2Al(n-n ) V 4n (n - n Q ) ' 



(46) 



(47) 



where 7 = ^2I p /Aq > is the Keldysh parameter [T] discriminating between adiabatic 
tunneling (7 < 1), non-adiabatic tunneling (7 ~ 1) [9], and multiphoton ionization (7 3> 1). 



Eq. (47) is further rewritten as 



Xn 



'1 + 7 2 

1-C 2 



(48) 



where 



C = 2 ^-ie(-l,l], 

n 



(49) 



corresponding to the range n > no- In this range, Xn is always larger than 1, thus Xn{k z ) > 1, 



cf. Eq. (46). Therefore, the ionization time U in Eq. (45) must be complex, i.e. U = Retj + 



i Im £j , and we get two equations for Re tj and Im i j 



sin(o; Re £, =F 0) cosh(o; Im tj) = Xn{k z 
cos(o> Re U =F Q) sinh(o; Im tj) =0. 



(50) 
(51) 
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The corresponding solutions for Imtj 7^ are 



7T 



u Re U = - ± + 2nN, 
culmti = arcoshx n (fc 2 ), 



(52) 
(53) 



where N e Z is chosen such that a; Ret, lies in the interval of a laser cylce, that is between 



Ti and Ti, cf. Eq. (38). The complex time t{ can be interpreted as the time of entering 



into the barrier, while its imaginary Imij and real Retj parts are the tunneling time and 



the time of exiting the barrier, respectively [H]. With Eqs. (44), (46), (52), (53), and 



sinhx = \/ cosh 2 x — 1, the corresponding action S±(k,ti) at k = k n (Eq. 



39) is 



S ± {k,t 



i)\k=k n 



7T COS 8 . . 

n I — ± p + 2txN — — + in\ arcosh Xn(k z 

2 Xn{k z " 



Xn(k z 



(54) 



Only its imaginary part 



Im 5(k, t 



t)\k=k n 



i 



1- 



Xn(k z y 



(55) 



does not depend on the sense (±) of circular polarization, thus we can omit the index ± of 
the imaginary part of the action. It means that the imaginary actions are equal for right 



and left circular polarizations. Using Eqs. (39), (44)-(46), the absolute value of the second 



derivative of the action S±(k, U) in Eq. (42) at k = k n 



dt 2 



S(k, t) 



1 



nuj 



(56) 



k=k n ,t=U V Xn(k z 

is also independent of the sense (±) of circular polarization. Then, the expression for the 



probability of the n-photon process (|42|) is rewritten as 
\F n ±(k,u 



I 2 - 



lm(y±(ti))\ k=kn -2n(arcoshxn(fc z )-A/l-l/Xn(M 2 ) 



27171-^/1 - 1/Xn(k z y 



(57) 



Thus, the dependence of the ionization rate on the sense of circular polarization is only due 
to the prefactor 



\<l>lm{v±{ti))\l = u 



-\{v ± (ti) 2 + 2I p )0 lm (v ± (t i )Y 2 



k — kn 



(58) 
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cf. Eq. (|9]). Using Eqs. ([6]), (52), and (53), the initial electron velocity 



v±(ti)\ k=kn = v x ±{ti)\ k=kn e x + v y ±(ti)\ k=kn e y + k z e z 



(59) 



with x- and y- components 



v x ±{ti)\ k=kn = k p cos9 - A (xn{k z ) cos 9 + i^/xn(k z ) 2 - 1 sini 
v y±( t i)\k=k„ = kpsin9 - A (x n (k z ) sin 9 ± i \/ Xn(k z ) 2 - 1 cos^ 



(60) 
(61) 



specifices the required momentum of the initial wavefunction, i.e. (pi m (v±(ti)). The amount 
of this momentum depends on the orbital, in particular it is different for p + and p- orbitals. 



Because of the saddle point equation (43), i.e. v±(ti) 2 \k=k n +2/ p = 0, the initial wavefunction 
ipim{v±(ti)) must have the pole at \±(ti) 2 \k=k n = — 2/ p = — k 2 that yields non-zero prefactor 



\4 , im( v ±(ti))\ k=kn , cf. Eq. (58). For short-range potentials, it corresponds to the wavefunction 
in coordinate representation asymptotically far from the core [12], i.e. 



(flm{r) = C k ik 3/2 Y lm (9 r , r ), 

nr 



(62) 



with the constant C K i, depending oiik= \/2I p and / as well as details of the potential near 
the core. Using spherical harmonics 



Y, 



'21 + 1(1 -m) 
47r (I + m) 



PI 71 (cos 9 r )e im ^, 



(63) 



Fourier transformation pi), and 



iv±(ti)r = — iv±(ti)r (sin 9 v ±(ti) sin 9 r cos(4> v ±(ti 



+ cos 9 V ± (ti) cos 9 r ) (64) 



in spherical coordinates, we evaluate two integrals over r and 9 r with the help of the Bessel 
function and Ref. [29] to yield the intermediate result (see Appendix 4) 



^m(v±(ti)) = C Kl J '-ll—YUe^iUl^iU)) e~ u ^ 2 I dr^e- Kr J l+ll2 (v ± (t t )r). (65) 
V W*) Jo 

Expanding the Bessel function in Taylor series and using the Gamma function T(z), the 
integration over r is easily carried out. The resulting series is then compacted as the Gaussian 
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hypergeometric series (see Appendix 4) 



<Pzm(v±(*i)) 



C K i ( v±(ti) 



Y lm (8v±(ti),<f>v±(ti))e- il ^ 2 



(66) 



T(/ + 2) // Z_ 3 3 t; ± (tj) 2 
r(Z + 3/2) 2 X V 2 2 2 2 ' k2 



The hypergeometric series does not converge at the saddle point v±(ti) 2 \k=k„ = —k 2 , thus 
this series has the pole as expected above. Multiplying Eq. (66) by v±{t.i) 2 + k 2 yields the 
series (see Appendix 4) 



(V±(ti) + K )<Plm(v±(ti)) = C Kt 



2k ( v±(ti)\ l 



7T \ K 

A r(z + 2) . 

2 Z +! T(/ + 3/2) 2 1 V2 2'2' 2 



^M*i)>&±(*i))e _<W2 
I 11.3 v±(ti) 2 



(67) 



which is convergent at the saddle point, i.e. 



A r(Z + 2) /I I 1,3 . 
2W T(Z + 3/2) 21 ^2'2 2' 2' 



(68) 



Thus, the prefactor (58) is simplified as 



V ± (t 



|C^|%/2/ 



2vr 



" |^m(^±(^),^±(^))lLfc„ 



(69) 



cf. Ref. [12J. With Eqs. (57) and (63), we finally obtain the general result for the probability 
of the n-photon process for all atomic orbitals 

2 



|F n± (k,u;) 



2 \C Kl \ 2 ^2T p (2l + l) (l-\m\)\ 

k=K 167^1 - l/xn(k z ) 2 (I + \m\)\ 

- 2n farcosh Xn(k z )- y/l-l/Xn(k z ) 2 ) 



P, 



\m\ 



iky 



/2X 



im<f> v ±(ti) I 



(70) 



where cos 9 v ±(ti)\k=k n = v z /v±(ti)\k=k n = ±ik z / \f2I p was used. The square of the associated 



Legendre polynomials in Eq. (70) are equal to 1 for s orbitals, k 2 /{2I p ) for p Q orbitals, 
{k 2 + 2Ip)/(2I p ) for p± orbitals, and so on. The ionizaztion rates for orbitals with m = 
(e.g. s and po orbitals) are independent of the sense of circular polarization. For m ^ (e.g. 
p± orbitals), they depend on the polarization sense, solely due to the factor \e tm ^ v± ^\ 2 k=kn . 
For m 7^ 0, this factor is not equal to unity, because the so-called tunneling momentum 
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angle 4> v ±(ti), which is related by 



COS (j) v ±(ti 

sin 4> v ±(ti 



Vx±(tj) 

v P ±(U) 
Vy±(ti) 



(71) 
(72) 



is complex. With v p ±{U 



-(^ 2 + 2/ p ),7= J2T P /A , Eqs.(31), (46), (60), and (61) 



the factor \e im ^ v± ^\l =kn for m = ±1 is (see Appendix 5) 



D im(f> v ±(ti) I 



\cos (f) v ±{ti) + isgn(m) sin<f) v ±(t 



(73) 



2 X n(k z ) 2 (l T sgn(m)Vl-l/xA) 2 ) - (1 + 7>M>] * 



, (74) 



27 2 Xn(^) 2 (^ 2 + 2/ p ) 

Now, we can see that the ionization rates for a given circular polarization are different 
for orbitals with opposite quantum numbers m = ±1 (e.g. p± orbitals). Because of the 



term =)=sgn(m) in Eq. (74), the ionization rate for the p + (or pJ) orbital and right ciruclar 
polarization is the same as the ionization rate for the p_ (or p + ) orbital and left circular 
polarization, supporting the fundamental symmetry in electrodynamics. 



E. Accurate formulas for the time-averaged ionization rates for s and p orbitals 



Since the function \F n ±(k,u)\ k=kn (Eq. (|70|)) depends only on k z , the two integrations 



over k p and 9 in the formula for the time-averaged ionization rate (Eqs. (35) and (36)) are 
easily carried out and the result is simplified to 



'0 



w±(£ ,u) = 8n 2 / dk z \F n ± (k, u 



n>n,Q 



I'- 



(75) 



With Eqs. (70) and (74), the accurate formulas for the time-averaged ionization rates are 

^-2n(arcoshxn(fc 2 )--\/ 1 ~ 1 /Xn(fc z ) 2 ) 



w s (£,u) 



\C K o\ 2 



V \ ^ 

7T ^ 



OO ., i>U 



n>riQ 



n 



dk y 



^1 - l/Xn{Kf 



(76) 
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for s orbitals, 



w P0 (S,u) 



00 1 rkn r,2 -2n(arcoshxn(fc z )--v/l-l/Xn(fc z ) 2 ) 

= / ^ fc * e 



3 \CkI\ 



(77) 



for p orbitals, and 



3|C k i| 2 a/2/ p ^-^ 1 



167T7 2 14 — ' n / n 

n>no 



fc n -2n(arcosh Xn (fc z )--\/l-l/xn(fc z ) 2 ) 



[2 Xn (k z ) 2 (l =F sgn(m)v/l-l/x„(A; 2 ) 2 ) - (1 + 7 2 )n/n 



(78) 



for p± orbitals. 



F. Approximate formulas for the time-averaged ionization rates for s and p orbitals 



In Eqs. (76)-(78), the exponential function has the maximum at k z = 0, confirming our 



expectation that the electron leaves mostly in the polarization plane of the laser field, i.e. 
x/y-plane, and that the electron ionization along the propagation axis (z-axis) is suppressed. 
Therefore, we use Taylor series of the exponent at k z « up to second order, i.e. 

2 



1 f k \ 

arcosh Xn{k z ) - y/l - l/x n (k z ) 2 « arcosh^n - V 1 ~ l /xl + ^ l - Vxl ( -jf) ' 



(79) 



The /^-dependent prefactors in Eqs. (76)-(78) are then replaced by the non- vanishing lowest- 



order terms of the corresponding Taylor series at k z f» 0, i.e. 



y/l ~ 1/XniKY y/l ~ I/Xl 



(80) 



for s orbitals, 



y/l ~ l/Xn{k z ) 2 X/I^l/Xl 



(81) 



for po orbitals, and 



2 Xn {k z ) 2 1 =F sgn(m) y/l - l/ X n{k z ) 2 - (1 + 7>M) 



Xn{k z ) 2 y/1 ~ l/Xn{kzY 
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(82) 



^Xl (l T sgn(m)0 - Vxf) - (1 + 7 2 W™o 



xW 1 - VXr 



for p-t orbitals. For po orbitals, the k z - dependent prefactor in Eq. (81) has no zeroth-order 
term due to the existence of the factor k\. That means that the ionization rate for p 
orbitals in the polarization plane is zero because of the destructive interference coming from 
two phase-opposite lobes. It also concludes that the approximation k z « for p orbitals 
may be not very appropriate, i.e. the electron from the po orbital will leave with non-zero 
final momentum component k z ^ 0, i.e. parallel to the z-axis, due to its orbital shape. With 



the Taylor approximations (79)-(82), we evaluate the remaining integrals in Eqs. (76)-(78) 



over k y as 



dk z e a ™(fcn) 



a/7t k n erf (a r 
2a n 



(83) 



and 



2 I k z 



V / 7rerf(a n ) e a " 



2al 



(84) 



where a r 



^1 — 1/Xn) 1 ^ 4 - For u) <C I p , the minimal number of photons Uq (Eq. (31)) 



must be very large. It follows that n > n ^> 1, hence Xn ^ 1 (cf. Eqs. (48) and (49)) 



and a n 3> 1. In the limit a n — > oo, the error function erf (a n ) and the exponential function 



e a ™ tend to unity and zero, respectively. For oo <C I p , we use Eqs. (76)-(84) and Eqs. (31), 
(46), A = y/2Ip/y, i.e. k n = nwy/(xn a/2/ p ) and uj = I p (l + 7 2 )/(n 7 2 ), to obtain the 



approximate formuals for the time-averaged ionization rates 



w s (S, oS) 



|C k q|%(1 + 7 2 
4 v /7rn 7 



E 

n>no 



-In (arcosh Xn~ \A~VXn) 
nXn(l ~ 1/Xn) 3/4 



(85) 



for s orbitals, 



w po (£,co) 



3|q 1 |%(l + 7 2 ) 3 

3-3 



32^y^rnQ'y 



E 

n>riQ 



n e 



- 2n (arcosh \ n - i/l-l/Xn) 



xl{i-i/xl? /4 



(86) 
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for po orbitals, and 



3|C Kl |%(l + 7 2 



7r n 7 



n>no 



Xn 6 



-2n(arcosh Xn-y/ l-l/xl) 



(87) 



[a/1 - 1/Xn T (2n /ri - 1) sgn(m) 



for p± orbitals. These expressions Eqs. (85)-(87) depend on Xn and n. If we use Eqs. (48) 



and (49) as well as arcoshx n = artanh-y/l — l/x 2 > then we obtain alternative expressions 
for the ionization rates depending on (, i.e. 



\r< 12 J / 1 \ !/ 2 

w (S,u) = fj- 1 + 

D'H)VR(^) 

n>n VS 1 ' / 



i±7iy /4 P -^H h \/s?-^ 

C 2 + 7 2 



(88) 



for s orbitals, where it coincides exactly with Eqs. (68) and (69) of Ref. [12] . 

3/2 



w Po (£,u;) 



3|C re i| I p f 1 
16\/2^no /2 V 7 2 



n>n 



1 7 2 ^ 5/4 e -^(a rt an h ^-^ 



for orbitals, and 



(89) 



«4 (£> W ) = I ^ 3/2 1 1 



8V27T n 



3/2 oo 



C 2 + 7 2 



7^ 



=F C sgn(m) 



i /i+ 7 2 v /4 



VRU 2 + 7 



rR ^ artanh ^ i-±^- - J ^-±^- 



(90) 



for p± orbitals. Since \/ (C 2 + 7 2 ) /(I + 7 2 ) > |C|, we recognize in Eq. (90) that the n- 
photon ionization rate is maximal for ^psgn(m) = 1 if ( > and for ±sgn(m) = 1 if ( < 0. 
Therefore, for ( > (low photon and kinetic energies) and e.g. for right circular polarization, 
the rate for p_ orbitals is larger than the one for p + orbitals. For ( < (high photon and 
kinetic energies), however, the rate for p + orbitals is larger than the one for p_ orbitals. For 
( = 0, corresponding to the photon energy nu = 2uquj = 4U P + 2I p and electron kinetic 
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energy k\j2 = (n — no)u> = hqu = 2U P + I p , the ionization rates for both p± orbitals are 



identical. By the way, we would like to stress that the ionization rate for p orbitals (Eq. (89)) 
is very small compared to the rates for s and p± orbitals, due to the ionization supression 
in the polarization plane. 



G. Simple formulas for the time-averaged ionization rates for s and p orbitals 

To obtain the simple analytical expressions for the ionization rates, the summation over 



n-photon processes in Eqs. (88)-(90) can be replaced with integration over C, i.e. 



£ « / dn = 2n — ^— . (91) 
n>no Jn ° J - 1 V ^ 

For oj <C I p , i.e. n 3> 1, the saddle point method for integration over ( is then applied, 



where the exponent in Eqs. (88)-(90) 



has a unique maximum at C = Co (7)- This maximum is determined by the saddle point 
equation 

d 



= 0, (93) 

C=Co 



that yields the transcendental equation for Co (7) 



or equivalently 



artmlM /^ = r^Jp^, (94) 

W 1 + 7 2 1 -CoV 1 + 7 2 



artanh^l - l/ X l_ = \-^~^l-l/xl^ (95) 

" * "-max lb v 



where Xn max = a/(1 + 7 2 )/(l - Co) ( cf - E Q- (j48|)) and n max = 2n /(l + Co) is the number of 
photons for which the n-photon ionization rate is maximal, corresponding to the electron 
kinetic energy k\^j2 = (n max — n )u = (2U P + I p )(l — Co)/(l + Co)- Since the solution of 
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Eq. (94) is in the positive range Co (7) £ [0, 1], the maximum of the (n max -photon) ionization 
rate for right (left) circular polarization is dominated by the electron ionization from the p_ 
(p + ) orbital. In the adiabatic limit (7 <C 1), the solution is approximated as Co(7) ~ 7 2 /3 
(see Ref. [12] and Appendix 6), corresponding to the electron kinetic energy fc 2 max /2 rj 
(2U P + I p )(l — 27 2 /3) ~ 2U P + Ip/3, whereas in the non-adiabatic limit (7 ^> 1), it yields 



Co (7) ~ 1 — 1/ In 7, see Ref. [12J. The exponent (Eq. (92)) at the saddle point Co (7) is 

2£ 



S(Co,7) 



3S 



(96) 



where S = (2/ p ) 3 / 2 and 



0(7) 



3Co 



7 2 (1-C 2 ) 



v/(l+7 2 )(Co 2 /7 2 + l) 



(97) 



and it does not depend on orbitals. We also need the second derivative of the exponent 
S"(Co,7),i-e. 



2 S((,j) 



C=Co 



4n (2C 2 + Co 2 7 2 + 7 2 ) Co 2 + 7 2 
;i + Co) 3 (l-Co) 2 (Co 2 + 7 2 )V !+7 2 



(98) 



to apply the saddle point method according to 



2vr 



i ^/(C)e^^/(Co)e^"> v _ s „ (Coi7) 



(99) 



Therefore, with these equations and Eqs. (88)-(90), we obtain the compact expressions for 
the time-averaged ionization rates 



w s (£, u) = \C K0 \ 2 I P £r ^(7) e"^ 9(7) 
Icq 



(100) 



for s orbitals, 



w P0 (£,cu) 



for p orbitals, and 



|C«i|%^^°(7)e-^ 9{7) 



|ai| 2 / P ^^(7) e -- 9(7) 



(101) 



(102) 
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for p± orbitals, where 




h,w = (1 - &) Wit)S^a + i) (103) 



W(7) = k'( 7 )^ d-cJ^j (104) 



^'^w^S W^^^l • (105) 



For s orbitals, Eqs. (97) and (103) coincide exactly with Eqs. (73) and (74) of Ref. [12]. In 
the adiabatic limit 7 < 1, i.e. Co(7 < 1) ~ 7 2 / 3 (l - 287 2 /45) (see Ref. [V> and Appendix 
6), the exponent and prefactors in Taylor series up to second order in 7 are 

3(7 < 1) « l- 7 2 /15 (106) 



^(7<l)«#-fl + ^) (108 



/i s (7 < 1) ~ 1 (107) 

2£ V 9 

3 -Y 2 
^( 7 < 1) « -T7Sgn(m) + ^-. (109) 

Therefore, the ratio of the ionization rates for p± orbitals 

t + t" (7<<1) (110) 

is always larger than 1 for right circular polarization and smaller than 1 for left circular 
polarization. That means, for e.g. right circular polarization, the ionization from p_ orbitals 
is more preferred than the ionization from p + orbitals. In the adiabatic case 7 = 0, i.e. 
the tunneling ist much faster than the rotation of the laser field, the ionization rates for p± 
orbitals are equal as expected. In the non-adiabatic limit 7 ^> 1, i.e. Co (7) ~ 1 ~~ V l n 7> the 
exponent and prefactors are approximated as 

3 In 7 . 

7 »l «_ -L (ill) 
27 

ftS(7 11 " (to^Ti (U2) 

"•^"koZf (U3) 
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7 



IT 1 



4(ln 7 ) 1 /2 

Thus, the ratio of the ionization rates for p± orbitals is 



In 7 



1 2 



sgn(mj 



(21n 7 ) ±2 ( 7 »1), 



(114) 



i.e. for right circular polarization, the ionization rate for p_ orbitals is much larger than the 
one for p + orbitals. 



III. RESULTS AND DISCUSSION 

For application, we use Kr atom in the ground state with ionization potential I p = 0.5 a.u. 
An infrared circularly polarized strong laser field with typical experimental parameters of 
the laser frequency lo = 0.057 a.u. (800 nm) and the laser amplitude £ = 0.06 a.u. (/ = 
2.5 • 10 14 W/cm 2 ) ionizes an electron from the 4p valence orbital of the Kr atom. In this 
experimental example, the Keldysh parameter is 7 = oj/£^J 2I p = 0.95, thus an 4p electron 
tunnels the ionization barrier non-adiabatically with respect to the rotation of the electric 
field. Although this is indeed non-adiabatic tunneling, many previous theoretical works 
are based on adiabatic approximation that cannot predict the difference of the ionization 
rates for p + and p_ orbitals in circularly polarized laser fields. With our analytical formulas 
derived in Section II, which are beyond the original work (TSJ H5] , we present the results in 
Figs.[Ifj3] and show that the ionization rates for p + and p_ valence orbitals are indeed very 
different, supporting our physical interpretation in the previous work [TO] . We also show 
the results for the ionization rates for po orbitals in Figs.[TJ-[3j Although the results for the 
ionization rates for s orbitals have nothing to do with the ionization of the Kr atom, we 
would like to include these results in Figs.[TJ-[3]as well, but these results cannnot be compared 
with the ones for p orbitals due to the generally different factors C K \ for s and p orbitals 
depending on the model system. However, in all calculations, we have used C K \ = 1 for 
simplicity. 

In Fig.[TJ the results for ionization rates depending on the laser frequency up to u = 
0.12 a.u. for laser amplitude £ = 0.06 a.u. are shown. The orange, green, blue, and red curves 
correspond to the rates for s, p , p + , and p_ orbitals, and the associated solid, dashed, and 
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dotted curves correspond to the accurate (Eqs. (76)-(78)), approximate (Eqs. (88)-(90)), and 



simple (Eqs. (100)-(102)) results, respectively. For s orbitals, the approximate and simple 



results coincide with each other very well within grafical resolution. The accurate results 
for s orbitals are a little separated from approximate and simple results, mainly due to the 



integral approximations (83) and (84) for u <C I p . For po orbitals, the ionization rates are 



very small compared to the ones for p± orbitals, supporting our thoughts in Section II. E, i.e. 
the destructive interference coming from two phase-opposite lobes of the p^ oribtal causes 
ionization supression in the polarization plane perpendicular to the orbital nodal axis. Again 
by more precise inspection, the approximate and simple results are similar while the accurate 
results are a little separated from the approximate and simple results, again mainly due to 



approximations (83) and (84) for oj <C I p . While the ionization rates for s and po orbitals 



do not depend on the sense of circular polarization, the rates for p + and p- orbitals do. For 
right circular polarization, the rates for p_ orbitals are larger than the ones for p + orbitals by 
the ratio up to 6 for large frequencies. For left circular polarization, the physical behaviour 
is reversed according to the fundamental symmetry in electrodynamics, i.e. the ionization 
rates for p± orbitals and left circular polarization are equal to the ones for p T orbitals and 
right circular polarization. In the adiabatic limit 7 = u = 0, the rates for both p± orbitals 
are exactly equal as already predicted in many previous theoretical works based on adiabatic 
approximation. For p± orbitals there are some (but not large) deviations between accurate, 
approximate, and simple results in particular for large frequencies or equivalently for large 
7, but in the adiabatic (7 <C 1) and non-adiabatic (7 ~ 1) tunneling regimes, these three 
results for p + and for p_ converge well. The small differences for large frequencies are not 
only due to the saddle point method (which is only applicable for low frequencies) but also 



due to the integral approximations (83) and (84) for uj <C I I 



Fig. [2] shows the monotonically increased ionization rates for s and p orbitals versus laser 
intensity I = c 2 Eq£ 2 in the range from 2 • 10 13 W/cm 2 to 2 • 10 14 W/cm 2 (corresponding 
to the laser amplitude S in the range from 0.0169 a. u. to 0.0534 a.u.) for laser frequency 
u> = 0.057 a.u. (800 nm) in logarithmic scale. In this figure, the accurate, approximate, and 
simple results coincide within graphical resolution. As already explained above, the rates for 
Po orbitals are small compared to the rates for p± orbitals. For right circular polarization, 
the rates for p_ orbitals are larger than the ones for p + orbitals. The corresponding ratio is 
large (small) for low (high) laser intensities or equivalently for large (small) 7. 
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Fig. [3] shows the photoelectron spectra for s and p orbitals, right circular polarization, 
laser amplitude £ = 0.06 a.u., and laser frequency u = 0.057 a.u. (800 nm). The curves are 



calculated using time-averaged n-photon ionization rates w n+ (£,uj) (Eq. (36)) versus final 



electronic kinetic energy at the detector k 2 n /2 (Eq. <m}). Since there are no sums in simple 



results (Eqs. (100)-(102)), only the corresponding accurate (solid) and approximate (dahsed) 



spectra (cf. Eqs. (88)-(90) and Eqs. (100)-(102)) are presented in this figure, as well as the 
spectra for total p orbitals according to w^ + (£, u>) = w^ + (£, u) + w^.(£, u) + Wn+(£, oj). In 
fact, the spectra are different for electrons coming from different orbitals. For right circular 
polarization, the ionization from p_ orbitals is dominant, but there is a unique kinetic energy 
for which the n-photon ionization rates for p + and p_ orbitals are equal. This is the final 



kinetic energy 2U P 



1.05 a.u. for the approximate results. For the accurate results, 



the intersection of the spectra for p + and p- orbitals lies at the energy a little more than 
2U p + Ip. Below (above) this unique electronic final kinetic energy, the ionization rates for p_ 
orbitals are larger (smaller) than the ones for p + orbitals. Therefore, the final kinetic energy 
indicates the strength and the direction of the ring current [2HJ EZ] generated in the ion, 
measured in correlation with the electron. Low energy electrons correlate to the ions with 
positive ring currents, while higher energy electrons correlate to the ions with negative ring 
currents. Furthermore, the locations of the maxima for s, po, and total p orbitals are similar, 
whereas the ones for p^ orbitals are shifted to lower energy and the ones for p + orbitals are 
shifted to higher energy. The reason is that the counter-clockwise ("positive", right) sense of 
circular polarization drives the electron from p_ and p + orbitals with clockwise ( "negative" ) 
and counter-clockwise ("positive") azimuthal velocities, yielding smaller and larger kinetic 
energies, respectively. In the adiabatic limit 7 <C 1, all photoelectron spectra have its 
maxima at 2U P ~ 0.55 a.u. and in this case the spectra for p + and p_ orbitals are identical. 



IV. CONCLUSIONS 

We have validated the approximations used in our previous publication to derive simple 
formulas for ionization from different sub-states of the p orbitals. We extended the PPT 
theory to strong field non-adiabatic ioniziation for valence p orbitals and we derived the 
corresponding ionization rates in full analytical form. Due to the existence of the complex- 
valued tunneling angle in the prefactor of the ionization rate, the rates are different for 
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degenerate p + and p_ orbitals and depend on the sense of rotation of the circularly polarized 
laser fields. Strong field ionization preferentially removes a counter-rotating electron. As 
expected ionization rates for degenerate p + and p_ orbitals are significantly larger than the 
rates for po orbitals due its orbital symmetry. We have also demonstrated that ionization 
rates and electron spectra obtained in this approach are gauge- invariant, unlike the results 
of the strong field approximation. 

An important extension of this work is the consideration of the electron spin [30] to 
describe electronic ring currents in the ion, which couple electronic spin and orbital degree 
of freedom. Other possible extensions of this work include the theory of the non-adiabatic 
ionization for p orbitals in circularly or elliptically polarized laser fields and static magnetic 
fields, see also works for s orbitals [3U [32] . 
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APPENDIX 1 



Here, we evaluate the 0-integral in Eq. (26) 



2t 



Jj = / d0e i ^- fcl ") M ^ fa »- fc ^^«(/±(A;i p ,0 1 ,0,t)+/ ± (A; 2p ,0 2j 0,t)). (115) 



o 



Using the expression for the function f±(k p ,6,4>,t) (Eq. (13)) and the Euler's formula, the 



pref actor of the integrand in Eq. (115) is rewritten as 



f±(k lp ,e 1 ,<p,t) + f±(k 2p ,e 2 , ( f>,t) (lie) 

\ {ki P e^ + k 2p e ie > ± 2iA e ±iut ) e"* + \ {k Xp e~^ + k 2p e~ w * =F 2zA e^) 



hence 

»2tt 



1 f 

h = - (k lp e idl + k 2p e iB2 ± 2iA e ±iuJt ) / fy e ^ e WV^-ki*)™<H-V*y-ki y )dn4>) (117) 
1 /" 2l 



( e i<f> e ip{(k2x—kix) cos <p+(k2y-ki y ) sin< 



We use the substitution |k 2 || — ki||| sin<// = (A; 2z — k\ x )cos<f) + (k 2y — ki y )sxn<f), which is 



satisfied by the relations (sin = a/T — cos 



2 . 



(fc 2 x - ki x ) sin0' + (£; 2? , - h y ) cos0' 

cos0= — (118) 

l k 2|| - ki|| | 

oin^ - ^ 2y - kly ^ Sin ^ - ^ 2a ' - klx ^ COS ^ (119) 
1^211 - km I 



With Euler's formula, we obtain 



e ±i<f> _ Ti(k 2x - fcig) + (k 2y - kiy) ^ ±i(f) , {120) 
|k2|| - kiy | 

= T< **f -**f (121) 

k 2 || — kiy | 



and <i0 = c?0'. The integral I\ is then reexpressed as 



h = t {k lp e^ + k 2p e ie > ± 2iA„e ±iwt ) f lp ^ - [** d<j>' e -^>N r ki||l W ( i2 2) 

2 |k 2 || — kiy | J 

_1 (fc e - ifl i + fc 2p e- 492 =f 2zA e^) fc2 f ^ Z [** d<f / e H/ Jp^-^ ™<f>\ 

2 k 2 ||-ki|i J 
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With the definition of the Bessel function of the first kind 

i p2ir / i\n p2tt 

J n ( x } = _j d e -in<l> e ix sin </> = ^_)_ j d( p^n4> e ix S m<p^ (123) 

the integral is evaluated as 

h = ™ [{k lp e i9 i + k 2p e ie > ± 2iA e ±i " t ) (k 2p e' m - k lp e- idl ) (124) 
+ (h P e-^ + k 2p e^ T 2iAoe**) (k 2p e^ - k lp e^)] 

| K 2|| — | 

which is further simplified to 

h = 2m {k\ p - k\ p - 2k 2p A sm(ut T 2 ) + 2k lp A sm(ut T Jl ^ k ^~ k f\ (125) 

| K 2| ~~ I 

cf. Eq.pl. 



APPENDIX 2 



In this appendix we prove the relation (Eq. 28 ) 

lim J dk\\ ^(k,|) pJl( k pk ^ = 2tt J dk|| ^(k||)5(k||) (126) 



with the arbitrary function g(ky) and ku = k p as 



/n T Ink 1 f P 00 

dk||p(k||) : = lim / d6 I dk p g{k p cos 9 e x + k p sin 9 e y )pJi(pk p ) (127) 

k\\ p-+°°Jo Jo 

f 2n f°° fx X \ 

lim / d9 dx g - cos9e x -\ — sin6>e,, J\(x) (128) 

p-Wo Jo \P P ) 

»2lT 



PZ7T poo 

/ d9 dxg^J^x) (129) 
Jo Jo 

poo 

2irg(0) / dxJ x {x) (130) 
Jo 

2*9(0) (131) 
2tt j dk||<7(k||)<y(k||). (132) 
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APPENDIX 3 



In this section we will evaluate k + - and ^--integrals in Eq. (29). First with Eq. (31), we 



rearrange the denominator of the /c + -integrand in Eq. (29) as 
'1 



16 



(k + - k^y + ky + k; + AQ + 2I P - 2n^ + iS 
(k+ + /c_) 2 + kl + k 2 z + Al + 2I P - 2n 2 u - iS 



[(k+ - k^f - AKi + US] [(k + + k.f - 4K 2 - US] 



(133) 



where Ki j2 = 2(ni j2 — n )u — k 2 . — k 2 z . It has four zeros k + = ki±, k 2 ±, i.e. 



k 1± = ±2 x /K 1 -iS + k. 
k 2± = ±2JK 2 + iS-k. 



(134) 
(135) 



Using Taylor series lim^o \JK ± iS = lim^o (y/K ±iS/ (2\/K)j , the zeros are rewritten as 



h± = ±2\/Ki + k_T 
k 2 ± = ±2\/K~ 2 -k_ ± 



iS 
iS 

7Wi 



(136) 
(137) 



Thus, the fc + -integrand in Eq. (29) is 

"1 



(k+ - k-) 2 + ky + k 2 z + A 2 + 2I P - 2n x uj + iS 

(k + + k-Y + k 2 y + k\ + A\ + 2I P - 2n 2 u - iS 
16h±{k + ) 



(138) 



h±{k+) 



(k+ - ki+)(k+ - fei_)(fc+ - k 2+ )(k + - k 2 _) ' 
The corresponding integral is then evaluated using the residue method as 



(139) 



/oo 
dk + ~h±(k + ) = 2m 6(Im k) Res ~h±(k) 

■OO 7. . I. 



(140) 



k=ki±,k2± 
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where the residues are calculated as 



Kesh±(k) = lim {k + — k)h±(k + ), 

fc+— >k 



(141) 



in particular 

Res h±(ki- 



8h±(k 



1+) 



2^K[ 



is 



2JK X - 2WK 2 + 2k_ - 



id 



iS 



-1 



(142) 



Res h±{k\ 



2^/K x +2^K 2 + 2k- 
8/i±(*i-) 



id id 
+ 



2^K[- 



iS 

/Ki 



2JK l + 2JK 2 - 2k_ - 



id 



iS 



(143) 



Res h±{k 2 



2^K X - 2y K 2 - 2k _ 
8h ± {k 2+ ) 



i5 



i5 



2^K~ 2 



iS 

/Ki 



s/Ki VK 2 
2JK[-2JK~ 2 + 2k 



i5 



iS 



(144) 



/ — i — id id 
2^/K[ + 2^/K^ - 2k_ — + 



Res h±{k 2 -) 



8h±(k 2 -) 



2JK~ 2 + 



iS 

/Ki 



2JK X + 2JK 2 + 2k- - 



id 



+ 



id 

Wo 



(145) 



2WK X - 2JK 2 - 2k. 



iS 



id 



IK, 



For S = 0, k_ = and ^ n x ^ n 2 ^ 0, i.e. ^ K x + K 2 ^ 0, all four residues (142)-(145) 



are finite. Therefore, the /c_-integral in Eq. (29) would be zero due to the existence of the 



factor in the integrand. Thus, the condition n = n x = n 2 , i.e. K — K\ — K 2 , must be 



satisfied and the residues ( 142 )-( 145 ) are simplified to 

2h±{k u 



Res h±(ki + ) 
Res h±{k x _) 
Res h±{k 2+ ) 





fK- 


iS \ 

VkJ 


(2yfK + kSj 
2h±(h_) 


(*_ 


id \ 




fK- 


iS \ 

Vk) 


(2y/K - fc_) 
2h±(k 2+ ) 




+*) 



2v^+^) (2VK-kS) (k. 



iS 

sfK 



(146) 
(147) 
(148) 
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For K < 0, only two poles k\ + and k 2 + have positive imaginary parts, therefore only two 
corresponding residues contribute to the fc + -integral (140), but in the limit 5 = and 
k- = 0, the integrand of the ^.-integral in Eq. (29) or the sum of these two residues times 



is exactly zero. Therefore, we consider only the remaining case K > 0, where only two 



poles ki_ and k 2 + have positive imaginary parts. Multiplying Eq. (140) by applying 



lim^o, and using residues (147) and (148) yields 



lim k_ 

5^0 



dk + h±(k + ) = 2ni 



h±{-2VK + fc_) - h±{2^K - k_ 
K ( 2\fK - k_ 



(150) 



Inserting this result of the fc + -integral in Eq. (29) and carrying the simple integration over 



k- yields the desired result (n = n\ = n 2 > n for K > 0) 

h±{2yfK) - h±(-2^K) 



^ poo PC 

lim J± (p,t) = tc \ / dk y I 

r „^„. J — oo J — c 



dk? 



n>n " ~°° 
,2 jl2 _ ?„2 /„ 2 1,2 



K 



(151) 



where K = 2(n — n )u; - k 2 - k 2 = k^ - k 2 - k 2 z , cf. Eq. (32) 



APPENDIX 4 



Here, we derive the wavefunction in momentum representation <^z m (v±(tj)). Using Fourier 



transformation pi), wavefunction in coordinate representation (fi m (r) (62), spherical harmon- 



ics Y lm { 



(63), and equation for — iv±(tj)r (64), we obtain 



<Pim{v±{ti)) 



2t + i(i-m)\ r drre - Kr 

(27r)^y 4vr (/ + m)!i 

gin^e-^^^^Wco^ P/^cos^) 



2?r 



im<f>r p —iv± (U)r sin 0„± (ij ) sin 6> r cos(4> v ±(U) — <f> r ) 



(152) 



Using the substitution <fi' = <fi r — 4> V ±(U) — tt/2, we get cos(0„±(£;) — r ) = — sin0' and 
d(j) r = d<p' . With the help of the Bessel function of the first kind ( 123[ ), we evaluate the 
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^-integral as 

d6 r sin 6 r e~ rv± {h)r cos e "± {u) cos 9r P[ n (cos 6 r ) J m (v± (U)r sin 6 V± {U ) sin 6 r ) . 
Using Eqs. (19) and (23) of Ref. [29], the ^-integral is evaluated as 

d9 r S i n re -™±^>cos8 v± (u)cos8 r p™( cos r ) J m ( u± (^)r sin w± (t;) sin0 r ) 



2ir 

e-^-^PricosO^t^J^v^uy), (154) 



v±(U)r 



hence 



<Pim(y±(U)) =a il /^i| m (^(i ! ),^±(t ! ))e- fc/2 | drV^e-^Jt+va^i^Or), (155) 



cf. Eq. (|65|). With the Taylor expression of the Bessel function of the first kind 

{-If fx 



J «« = £ a „vi/,, x ^ ( 5 ) • < 156 > 



/3=0 

we have 



^ m (v ± (t 4 )) = C Kl] I^^Y lm (e v± (t i ),cj )v± (t i ))e- il ^ 2 

2^B\T(P + 1 + 3/2) { 2 J J, 



(157) 



-1V» A,,^\ 2 ^ +1 / 2 /-oo 



/J=0 

Using the definition of the Gamma function (z > 0) 

/>oo 

T(z) = / dzz*-^ - *, (158) 
Jo 

the r-integration is easily carried out 

rfrr 2 " +,+1 e-" = Q^±i±^, (159 ) 
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hence 



<£zm(v±(*i)) 



C K 



kI 



(v±(ti)\ l 



Y lm (Ov±(t i )Av±(t i ))e- u ^ 2 



£ 

0=0 



2k 3 V 2k y 

1 T(2/3 + / + 2) / ^±(ti) 2 



/3! r(/3 + Z + 3/2) V 4k 2 



Using the duplication formula for the Gamma function 

r(z)r(z + 1/2) 



r(2*) 



2 1 " 



■22 



7T 



for z = (3 + Z/2 + 1 yields 



<£im(v±(*i)) = C K/ 



2 /^(ii) 



7TK V K 



Yi m (9 v ±(ti),(j) v ±(ti))e 



-Utt/2 



£ 

0=0 



1 + Z/2 + l)T{(3 + Z/2 + 3/2) / v ± {U) 



0\ 



+ Z + 3/2) 



2\ P 



Using the Gaussian hypergeometric series 

r(c) 



2 F 1 (a,b; c; z) 



rpg + a)r((3 + b) 



E 



and Eq. (161) for z = Z/2 + 1, we obtain 



<&m(v±(ti)) 



/Vfc(t»)Y 



W \ 2k J 



Y lm (dv±(U),(j> v ±(U))e- U7T/2 



Eg+2) // Z_ 3 3 v±fo) 2 
T(Z + 3/2) 2 1 V2 2 + 2' + 2' k 2 /' 



cf. Eq. (66). Using Euler's hypergeometric transformation 



2 F 1 (a,b;c;z) = (1 - z) c a b 2 Fi(c - a, c - 6; c; z), 



the wavefunction is then rewritten as 



0Jm(v±(*i)) 



G 



2k 




1>±(tj) 2 + k 2 V vr \ k / 



«±(ti) V 



^m(^±(*i) > ^±(*i))e- <i ' r/2 



g T(Z + 2) /Z 1 I 3 t; ± (*i) 2 < 
2'+ 1 T(Z + 3/2) 2 1 \2 2' 2' 2' k 2 
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cf. Eq. (67). With the Gaussian theorem for hypergeometric series 

T(c)T(c-a-b) 



2 F 1 (a,b; c; 1) 



r(c-a)T(c-b) 



(167) 



and Eq. (161) for z = 1/2 + 1, the series in Eq. (166) is convergent for v±(ti) 2 = —k 2 , i.e 

2 l+1 T{1 + 3/2) 



2 F 1 



2'2 + 2' 2' 



7T T(l + 2) 



(168) 



cf. Eq. (68), therefore the wavefunction in momentum representation (166) has the pole at 

V^) 2 = -K 2 . 



APPENDIX 5 



Using Eqs. (|7l|), J72(, and v p ±{ti) 2 \ k=kn = -{k 2 z + 2I P ), the factor | e <"»*»± (*0 1 ^ for 



m = ±1 is rewritten as 



D im4> v ±(ti) I 



|cos^±(tj) + isgn(m) sin </>„±(£i)| 
k 2 + 2I P ' 



(169) 
(170) 



With Eqs. (60), (61), k p = \/k 2 — k 2 , and \e im6 \ = 1, it becomes 

AaXn(k z ) (l =F sgn(m)v/l-l/ X „(A; z ) 2 ) - jW^W z 



imtj> v ±(ti) I 2 



/c 2 2Ip 



(171) 



Using Eq. (46), we obtain 



Jm4> v ±(ti) I 



A 2 X n{k z ) 2 (l =F sgn(m)v/l - l/xn(& 2 ) 2 ) - nu; 



A 2 Xn (k z ) 2 (k 2 + 2I p 



(172) 



Using 7 2 = 2/ p /A 2 and Eq. i.e. n w = A 2 /2 + J p = J p /7 2 + 7 p = 7 p (l + 7 2 )/7 2 , we 
obtain the result for m = ±1 



im4> v ±(ti) I 



2 X n(A:,) 2 (l =F sgn(m)v/l-l/x„(A; z ) 2 ) - (1 + 7 >/n 



2 7 2 Xn(^) 2 (fc 2 + 2/ P ) 



(173) 



cf. Eq. ((74). 
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APPENDIX 6 



To obtain the simple analytic solution Co (7) of Eq. (94) 



l+l' 1 -Co V ! + 7 
in the adiabatic limit 7 <g 1, we use power series 



Co(7) = 5>^> ( 175 ) 



i=Q 

yielding 



1+7 1 - Li=o c i7 V 1+7 



In the Taylor series of Eq. (176) at 7 = 0, each coefficient of powers 7* must be zero. For 
zeroth order, we get 

artanhco — — = 0, (177) 

1 - c 

with the solution Co = 0. The coefficient for the first order is then automatically zero. For 
the second order, we get 

^4^4=0, (178) 
V 1 + C I 

thus c\ — 0. For other orders, we get 

\ ~ c 2 = (179) 
c 3 = (180) 

+ 

135 



" 8 +c 4 = 0, (181) 



thus C2 = 1/3, C3 = 0, and C4 = —28/135. The solution of Eq. (94) is therefore 



1- / \ 7 2 A 28 2 236 4 52i2 6 1 2 570692 8 , 

Co (7) = I 1 7 + 7 7 + 7 (182) 

sv ,; 3 V 45 ' 567 ' 18225 63149625 ' 
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and for small 7 1 



Co(7) ~ t- ( 183 ) 
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FIGURE CAPTIONS 



FIG. 1: Time- averaged ionization rates w + (£,u>) for s (orange), p (green), p + (blue), p_ 
(red) orbitals and right circular polarization versus laser frequency u for I p = 0.5 a.u. and 
laser amplitude £ = 0.06 a.u. The solid, dashed, and dotted curves correspond to the accu- 
rate (Eqs. (|76(-(|78[), approximate (Eqs. g-Q), and simple (Eqs. (fl00|-(fl02|)) results, 



respectively. Note that \C, 



kI 



1 are used. 



FIG. 2: Time-averaged ionization rates w + (£,co) for s (orange), p (green), p + (blue), p_ 
(red) orbitals and right circular polarization versus laser intensity / = c 2 Sq£ 2 for I p = 0.5 a.u. 
and laser frequency u = 0.057 a.u. (800 nm) in logarithmic scale. The solid curves corre- 



sponds to the accurate results (Eqs. (76)-(78|). Because of the logarithmic scale, the ac 



curate (Eqs. (76)-(78)), approximate (Eqs. (88)-(90)), and simple (Eqs. (100)-(102)) results 
coincide within graphical resolution. Note that \C K i\ 2 = 1 are used. 



FIG. 3: Time-averaged n-photon ionization rates w n+ (£,u) (Eq. (36)) or equivalently pho- 



toelectron energy distribution at the detector for s (orange), po (green), p + (blue), p- (red), 
and total p (brown) orbitals and right circular polarization versus final electronic kinetic 



energy k n /2 (Eq. (37)) for I p = 0.5 a.u., laser amplitude £ = 0.06 a.u., and laser frequency 



0.057 a.u. (800 nm). The solid and dashed curves correspond to the accurate (Eqs. ((76 



(78)) and approximate (Eqs. (88)-(90)) results, respectively. The spectra for total p orbitals 
are calculated according to u^ + (£,Co>) = w^ + (£,u) + w^(£,co) + Wn+(£,uj). Note that 



1 are used. The approximate results of the ionization rates for p + and p- orbitals 



are equal at the final kinetic energy 2U P + I p ~ 1.05 a.u, see text for discussion. In the 
adiabatic limit 7 1, all photoelectron distributions are peaked at 2U P ~ 0.55 a.u. and are 
the same for p + and p_ orbitals. 
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